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Collision process between breather and moving kink-soliton is investigated both analytically and 
numerically in Fermi-Pasta-Ulam (FPU) chains. As it is shown by both analytical and numerical 
consideration low amplitude breathers and soft kinks retain their shapes after interaction. Low 
amplitude breather only changes the location after collision and remains static. As the numerical 
simulations show, the shift of its position is proportional to the stiffness of the kink-soliton, what is 
in accordance with the analytical predictions made in this paper. The numerical experiments are 
also carried out for large amplitude breathers and some interesting effects are observed: The odd 
parity large amplitude breather does not change position when colliding with widely separated soft 
kink-antikink pair, whilst in the case of closely placed kink-antikink pair the breather transforms 
into the moving one. Therefore it is suggested that the "harmless" objects similar to the kink- 
solitons in FPU chains could be used in order to displace or move the strongly localized structures 
in realistic physical systems. In particular, the analogies with quasi-one dimensional easy-plane type 
spin structures are discussed. 



PACS numbers: 63.20.Ry; 05.45.Yv; 63.20.Pw 

I. INTRODUCTION 

Chains of classical anharmonic oscillators can serve as 
models for more complex physical systems which under 
definite conditions could be treated as one dimensional 
objects, e.g., optical fibers, magnetic film waveguides, 
quasi-one dimensional spin systems, DNA, ionic crystals, 
etc. Modeling various physical processes one can directly 
see consequences using computer simulations and com- 
pare them with the established analytical schemes. In 
the present paper it is proposed to model different non- 
linear processes in chains of coupled oscillators making 
simultaneous interpretations and predictions concerning 
real physical systems. 

Already the simplest model, as the one-dimensional 
chain of equal-mass oscillators exhibits following non- 
trivial phenomena such as energy equipartitionEJIi ap- 
pearance of n various patternstl and localizations^ (ei- 
ther movmeQ or statia3,El), different regimes of chaotic 
dynamicsEl,cl, etc. Therefore these classical systems could 
serve as tools for better understanding of nonlinear phe- 
nomena in completely different (on the first sight) many- 
body systems. For instance, invariance under the sim- 
ple symmetry transformation u n — > u n + const (u n is 
a displacement ofji-th oscillator) relates Fermi-Pasta- 
Ulam (FPU) chainHl (interparticle forces are functions of 
only relative displacements) with a wide class of systems 
with continuous symmetriesEij, e.g. quasi-one dimea,- 
sional easy plane ferromagnets-.and antiferromagnetaI3, 
ferrimagnetic spiral structuresE3 and even quantum Hall 
double layer (pseudo) ferromagnetstil. Such systems 
are characterized by the infinitely degenerated energy 
ground state. Spontaneous breakdown of the symme- 
try (by choosing a definite ground state) leads to the 
appearance of gapless Goldstone mode forming the kink- 
solitons in low energy limit. These localized solutions 



are well known for easy plane magnetic structures! 




and demonstrate similar properties as in the case of FPU 
chains!^ Jil 

The main difference between Goldstone mode kink- 
solitons and ordinary kinks of the models related to 
sine-Gordon equation (particularly, its discrete analogy 
- Frenkel-Kantorova modeIi3,EE) is that the formers do 
not carry topological charge. Besides-ithat, the kinks 
are believed to be the exact solutionsEia in FPU chain. 
Because of these circumstances it is expected that they 
should not decay itself and do not destruct other localiza- 
tions during the scattering process as long as no energy 
redistribution is required. In this connection it should be 
mentioned that FPU chain, the linear spectrum of which 
is bounded from above, exhibits another nontrivial solu- 
tion in high energy limit. This solution. represents intrin- 
sic localized mode (discrete breathcr)Ell,E3, which in a low 
amplitude limit could be conskLeced as particular case of 
semi-discrete envelope solitonnJ,a. Let us note a direct 
analogy of the above with quasi one-dimensional mag- 
netic systems w-hepe similar localizations have recently 
been discoveredE3j23 or predictedEZl. 

As it follows from the analytical and numerical con- 
siderations made in the present paper the kink-solitons 
are indeed "harmless": after interaction the shapes of 
both kink and breather remain unchanged. The collision 
only causes the shift of the position of spatially local- 
ized breather or its transformation into the slowly mov- 
ing one. In this connection let us make a comparison 
with strongly inelastic scattering procesa-between kinks 
and breathers of sine-Gordon equationc3j2j. Althogh, it 
should be mentioned that in the latter case the nonlin- 
ear objects are solutions of continuous equation unlike 
discrete FPU model considered in the present paper. 

For analytical consideration in weakly nonlinear limit 
the multiple scale analysis will be used in order to present 
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the quantitative picture for kink-breather collision. But 
firstly, well known solutions for kink and breather will 
be briefly redflriyed in order to introduce the method of 
calculations 



II. ANALYTICAL SOLUTIONS FOR 
KINK-SOLITON AND BREATHER IN WEAKLY 
NONLINEAR LIMIT 

The equations of motion of FPU oscillator chain are 
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where dots over u n express the time derivatives. Di- 
mensionless units are used so that the masses, the lin- 
ear and nonlinear force constants and the lattice spac- 
ing are taken equal to unity. The real displacements are 
expressed from dimensionless ones (u n ) by dividing the 
latters on the coefficient y K^/m where m is a mass of 
particle and K4 is a coefficient before the anharmonic 
quartic term. Thus if the nonlinear interaction is strong 
enough it is admissible to have large values of u n (e.g. 
!i n > 1) and this does not cause the scattering of neigh- 
boring particles. 

Firstly let us derive the kink-soliton solution by as- 
suming that u n smoothly varies in space-time. Then it 
is appropriate to introduce slow variables 



£i = e(n - vxt); 



n 



eH 



and denote 



u n = Pl(fl,Tl), 



(2) 



(3) 



where e is a formal small parameter indicating smallness 
or slowness of the variables before which it appears. Sub- 
stituting (||) into the motion equation (Q) and collecting 
the terms with the same order of e it becomes possible to 
treat the problem perturbatively. In particular, the ve- 
locity v\ is determined in the second approximation over 
e: 



±1. 



(4) 



Without restricting of generality let us further consider 
the solution with negative velocity v\ = — 1. Other solu- 
tion will be recovered simply changing the axis direction. 
Finally, in the forth approximation over e the following 
nonlinear equation is obtained: 
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which is-exactly integrable modified Korteweg-deVries 
equatiorJ23 for the function The equation (||) 

was derived for FPU chain in Refs.t3,E-D and finally leads 
to the kink- like solution for u n : 



which has a similar form as the kinks for the sine-Gordon 
equation but note that although the tails of the kink solu- 
tion (0) corresponds to the different ground states (u„ — 
and u n — 7r/v6 for n — > — oo and n — > oo, respec- 
tively) these ground states carry the same energy because 
of the mentioned symmetry u n — > u n + const. These 
kinks do not carry topological charge and as far as they 
connect degenerate ground states they could be called 
Goldstone mode kinks, ft should be also mentioned that 
in terms of relative displacements v n — u n+ \ — u n this 
object is discretized version of Korteweg-deVries soli- 
ton and therefore the definition of kink-soliton is usually 
used in literature for its identification. The similar local- 
ized objects could be created in magnetic structures with 
easy plane anisotropy where their appearance also is con- 
nected with the broken symmetry Goldstone mode. The 
transverse component of such magnetic localization (in- 
plane component) has a kink like form, while out of easy 
plane component represents ordinary Korteweg-deVries 
solitonfcajii 

The solution (||) is valid if one can neglect the higher 
derivatives. This could be achieved if the following con- 
dition is satisfied for the kink stiffness: 



6A 2 < 1. 



(7) 



Afterwards let us rederive the breather solution using 
multiple scale analysis presenting multiplication of 

harmonic oscillation and smooth envelope function: 



ln = -^ 2 (&,T 2 y +C.C. 



(8) 



where "c.c." denotes complex conjugation and new slow 
variables are defined as follows: 



£2 = e(n - vit); r 2 = e 2 t. 



(9) 



As far as only small displacements are considered it is 
natural to neglect the higher harmonics working in a ro- 
tating wave approximation. Carrying out the procedure 
similar to the previous case (collecting terms with the 
same harmonics and order of e) in the first order over 
e well known dispersion relation for linear excitations in 
FPU chain is obtained: 



= u>k = \/2(l — cosk). 



(10) 



In the second approximation the expression for group 
velocity is derived: 



i'2 



sink _ dujk 
dk : 



(11) 



and finally we get the nonlinear Schrodinger equation for 
the envelope function if2 in the third approximation over 
e: 



"" = ^1 = V 3 ( arctan 
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which permits bright soliton solution. Thus in terms of 
u n the envelope soliton solution (mcpjng with a group 
velocity u 2 ) is rederived (see e.g. RefH): 



Bcos(nk — ujkt) 



v 2 t)] ' 



(13) 



LOk — <^h 



(l + (3/16)^5- 



B< 1. 



The breather solution is obtained by setting v 2 = 0, 
therefore carrier wave number k = n (thus u> = 2) 
should be considered according to the relations ( [Toj ) and 
(pd|). Thus we get the expression for the low amplitude 
breather solution: 



Bcos(nn-2t- (3/2)B 2 t) 



B < 1 (14) 



which coincides with the corresponding breather solution 
obtained in Ref.t 



III. INTERACTION BETWEEN 
KINK-SOLITON AND BREATHER 

A. Analytical Results in Weakly Nonlinear Limit 

Now let us start the main task of the paper: analytical 
description of kink-breather interaction. Keeping in mind 
that in absence of either kink or breather one should come 
to the solutions or (^|), respectively, I am seeking 
for the solution in the following form (using again the 
rotating wave approximation): 

. (15) 

where the following choice for slow space-time variables 
is made: 



6 =e{n + t)-e 2 * 1 {&,T 2 ), 



n = e% 



£ 2 = erc-£ 2 # 2 (£i,n), r 2 =e 2 t. 



(16) 



Here the phase and argument shifts are introduced in 
order to decouple nonlinear equations. Substituting jl5| ) 
into the initial equation of motion for FPU chain (Q) wc 
get in the forth order over e for zero harmonic and in the 
third order over e for the first harmonic the following two 
nonlinear equations: 
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The variables £i and £ 2 in Exp. (h6h are chosen such 
that group velocities for noninteracting kink-soliton and 
breather (with carrier wave number equal to n) are vi = 
— 1 and v 2 = 0, what itself guarantees the satisfaction of 
the motion equation (Q) in the lower orders over e. 
By letting 
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we come again to the equations (||) and (12) for kink- 
soliton ipi and breather ip 2 (with carrier wave number 
k = 7r) . The choice (|l9|) physically means that the inter- 
action effects reduce only to the phase shifts of solitons 
while the solitons' profiles remain unchanged in the lead- 
ing approximation. 

Finally, in the fourth approximation over e for the first 
harmonic the following equality is derived: 



3 / dcpx 
2 UCi 



(20) 



According to the last relation breather acquires group 
velocity during the interaction process, but as the kink 
passes it stops. The shift of breather's position could be 
calculated from the following simple relations: 



h = 



v 2 dt — 



dt 



dt = 



.( 21 ) 

thus shift is always positive, i.e. breather will be shifted 
oppositely to the kink's propagation direction (let us re- 
mind that group velocity of the kink has been chosen to 
be negative) irrespective of the sign of A, therefore the 
shift is the same for both kink and antikink case. 

Denoting by t\ and t 2 the times needed for kink to 
travel from the one side of the chain to another in pres- 
ence or absence of the breather on the way, one can cal- 
culate the difference At = t 2 — t\ using similar to ( |2l| ) 
relations. Thus one gets: 



(22) 



At = 2V6|B| 



The physical meaning of the expressions (|21j) and ( J22 
could be simply understood mentioning that in case of 
weakly nonlinear solitons' interaction the group veloci- 
ties of the solitons change only during the interaction 
process. Particularly, the breather acquires the nonzero 
velocity while interacting with kink-soliton. Simultane- 
ously, during the same small time period the velocity of 
kink-soliton becomes larger than in case of its free propa- 
gation. These circumstances cause the shift of low ampli- 
tude breather position and on the other hand the earlier 
arrival of the kink soliton at the left side of the chain (see 
Fig. 1). 
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teracton becomes again static. According to Fig. 1 the 
shift is equal to 1% s» 1 as it expected from formula j2l|), 
note only that as four localized objects with the same 
stiffness are used (two kink-antikink pairs) the value ob- 
tained from the expression ( ^l| ) should be multiplied by 
the factor 4. The interaction also causes acceleration 
of the kink-antikink pair (propagation velocity of kink- 
antikink pair is larger during the interaction process in 
comparison with free propagation) as it follows from the 
formula (122 
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FIG. 1: Shapes and locations of kink-antikink pairs and low 
amplitude breather (a) before and (b) after collision. Kink- 
antikink pairs move from right to the left. Insets show the 
enlarged view of the breather. u„ is the displacement of n-th 
site. 



B. Numerical Experiments 

As it was mentioned above, direction of the shift of 
the breather position does not depend on whether the 
kink or antikink participate in the collision. Therefore in 
order to increase the interaction effect two kink-antikink 
pairs are used for collision with breather. In the nu- 
merical experiment the soft kink solution (^) and low 
amplitude breather ( [l4"| ) are put apart from each other 
in FPU chain with pinned boundary conditions (see Fig. 
1). Numerical experiment fully confirms analytical pre- 
dictions: the nonlinear objects for which analitycal re- 
sults (^) and ([l4]) are valid [i.e. the conditions (Q) and 
( |l4| ) are satisfied], behave in full accordance with formu- 
las (|2l]) and (p2|). Particularly, as series of numerical 
experiments show, the shift of low amplitude breather 
position is proportional to the kink-soliton stiffness and 
does not depend on its own amplitude. In Fig.l the col- 
lision process between two kink-antikink pairs with the 
same stiffness A = 0.2 and the breather with amplitude 
B = 0.1 is expressed. It is clear that the localized ob- 
jects behave as expected: they retain their shapes and 
low amplitude breather changes its position and after in- 




FIG. 2: (a) collision of widely separated kink-antikink pair 
(stiffness A — 0.25) with strongly localized odd parity mode 
(amplitude B — 1). (b) collision of the same breather with 
closely placed kink-antikink pair. Arrows show the initial 
position of the kink-antikink pair. 

Obviously, the weakly nonlinear approach fails consid- 
ering large amplitude breather and/or stiff kink-solitons. 
Therefore it is expected that the results could be dif- 
ferent. Indeed, the following phenomena are monitored 
when colliding soft kink-antikink pair with strongly lo- 
calized odd parity mode. Breather with amplitude above 
B w 0.2 does not change the position if soft kink- 
solitons are used for collision. Moreover, breather does 
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not change the position if widely separated kink-solitons 
are used for collision. On the other hand, if there is 
a relatively short distance (by order of inverse stiffness) 
between kink and antikink the breather starts moving 
after collision and does not stop. This behaviour is ex- 
pressed in Fig. 2, where the strongly localazed odd parity 
mode (... 0, -1/2, 1, -1/2, ...) is placed in the center 
of FPU chain and different effects are monitored for dif- 
ferent separations between kink-antikink pairs (with stiff- 
ness A — 0.25). Apparently this some kind of "quantum" 
phenomenon is caused by the fact that single soft kink 
(or antikink) does not carry enough impulse to trans- 
form the static breather into the moving one, while the 
kink-antikink pair is able to do so. This statement is 
consistent, with the observation of chaotic breathers (see 
e.g. Ref.EJ) and with theJattice quantization procedure 
recently proposed in RefO. It should be mentioned that 
all localized objects retain their shapes after interaction 
and the acceleration of kink-antikink pair is still observed 
during the interaction process. 




FIG. 3: Collision of closely placed kink-antikink pair (stiffness 
A = 0.35) with larger amplitude (B=2) odd parity mode. 
Note that although the breather starts to move after collision 
(as in the previous figure) , further it is trapped by the lattice. 

Considering the breathers with larger amplitudes, nat- 
urally, the stiffer kink-antikink pairs will be required to 
displace the breather. A number of numerical experi- 
ments have been made for breathers with larger ampli- 
tudes. The general feature is represented in Fig. 3, where 
the collision of odd parity mode (..., 0, -1, 2, -1, 0, ...) 
and kink-antikink pair with stiffness A = 0.35 is demon- 
strated. Like the previous case the breather does not 
change the position colliding with widely separated kink- 
antikink pair. However, in contrast to the previous case, 
although the breather starts to move because of the col- 
lision with closely placed kink-antikink pair, after some 
time it is trapped by the lattice sites and stops. 

It should be mentioned that according to the analyt- 
ical results and numerical experiments the low ampli- 
tude breather acquires the group velocity only during the 
interaction process, while the large amplitude breather 



once starting to move does not stop (for intermediate 
amplitudes) or it will further trapped by the lattice sites 
(for larger amplitudes, approximately B > 1.4). 

In the present paper only the results concerning colli- 
sion process with odd parity mod* are presented because, 
as numerical simulations showEH, the even parity mode 
spontaneously starts to move and decays after finite time 
period, while odd parity mode remains well localized in 
the absence of collisions with other nonlinear objects. At 
the same time only soft kink-solitons are the subject of 
the study because stiff kinks sufficiently perturb back- 
ground and it is hard to see what causes displacement or 
pinning of the breather. 

Let us mention also that reflected kink-antikink pairs 
almost return the low amplitude breather to its initial 
position, while they are unable to change the picture in 
case of large amplitude odd parity mode. As it is seen 
form Fig. 2 moving odd parity mode with amplitude 
B=l does not react on the collision with reflected kink- 
antikink pair and remains moving with the same veloc- 
ity. On the other hand the reflected kink-antikink pair 
can not displace (see Fig. 3) larger amplitude breather 
(B=2). 



IV. CONCLUSIONS 

Summarizing, as it is shown above, kink-solitons could 
be used to displace or move the static localization with- 
out the destruction of the latter. The direct analogy 
with quasi-one dimensional easy plane magnetic struc- 
tures should be quoted again. The existence of contin- 
uous U(l) symmetry in magnetization vector space for 
this systems allows the presence of the broken symmetry 
Goldstone mode which because of the existing nonlinear- 
ity forms kink-soliton. These large wavelength nonlin- 
ear excitations have been studied for easy plane antifer- 
romagnets in presence of applied magnetic-field along 
anisotropic axisEj ancLin spiral structuresE3. On the 
other hand it is knowntHl that the band edge excitations 
in easy plane type antiferromagnets form stable intrin- 
sic localized spin wave modes (ILSM) having odd parity 
structure in large amplitude limit. Thus, a new theoret- 
ical study and corresponding realistic experiments could 
be planned in order to investigate and observe the effects 
caused by the interaction between magnetization kink- 
solitons and ILSM in quasi-one dimensional easy plane 
structures. As the analogy is almost straightforward it 
could be predicted that kink-solitons in the easy-plane 
magnetic structures should cause the same effect as in 
FPU chains. Particularly, they can displace or move the 
strongly localized objects without their destruction. 
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